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Abstract 

We show that quantum coherence of biomolecular excitons is maintained over 
exceedingly long times due to the constructive role of their non-Markovian 
protein-solvent environment. Using a numerically exact approach, we demon- 
strate that a slow quantum bath helps to sustain quantum entanglement of 
two pairs of Forster coupled excitons, in contrast to a Markovian environ- 
ment. We consider the crossover from a fast to a slow bath and from weak 
to strong dissipation and show that a slow bath can generate robust entan- 
glement. This persists to surprisingly high temperatures, even higher than 
the excitonic gap and is absent for a Markovian bath. 

Key words: Quantum entanglement in biomolecules, open quantum 
systems, Forster energy transfer in excitonic systems 



1. Introduction 

Quantum coherent dynamics at the initial stages of photosynthesis in 
complex biomolecular structures seems to promote the efficiency of energy 
transfer from the light-harvesting antenna complexes to the chemical reac- 
tion centers [H, 0, H 0]- This hypothesis has recently been boosted by ex- 
periments revealing long-lived quantum coherent excitonic dynamics in the 
energy transfer among bacteriochlorophyll complexes over a surprisingly long 
time of around 600 fs measured at 77 K In addition, electronic coherences 
between excited states in purple photosynthetic bacteria have been monitored 
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in a two-color photon echo experiment [2j. Both works lead to the conclu- 
sion that the collective long-range electrostatic response of the biomolecular 
protein environment to the electronic excitations should be responsible for 
the long-lived quantum coherence. Furthermore, the obtained time scales j^j 
for the short-time dynamics of the nuclear modes coupled to the excitonic 
states of two chromophores are almost identical. This points to the special 
and constructive role of the quantum environment for the photoexcitations. 
The often assumed coupling of the chromophores to fast and independent 
quantum baths does not hold in this case. In fact, the two chromophores are 
embedded in the same protein-solvent environment. These results corrob- 
orate experimental studies which show that energy transport sensitively 
depends on the spatial properties of the delocalized excited-state wave func- 
tions of the whole pigment-protein complex. In addition, there are reports of 
coherently controlled wave packet quantum dynamics artificially generated 
by laser pulses in the light-harvesting antenna of the bacteria Rhodopseu- 
domonas acidophila 

An appropriate theoretical description of the biomolecular quantum dy- 
namics has to account for the environmental time scales typically being of 
the same order of magnitude as the excitonic time scales [5j. This fact ren- 
ders the dynamics highly non-Markovian and rather elaborate techniques are 
required to address the entire cross-over from fast to slowly responding baths. 

Here, we perform a deterministic evaluation of real-time path integrals 
0, 0, S| where all non-Markovian effects are exactly included. We provide 
numerically exact results for the quantum coherent dynamics of photoexci- 
tations in coupled chromophores, where the time evolution of the protein- 
solvent bath happens on time scales comparable to the exciton dynamics. We 
show that the non-Markovian effects help to sustain quantum coherence over 
rather long times. Furthermore, quantum entanglement of two chromophore 
pairs is shown to be more stable under the influence of a non-Markovian 
bath. Even at high temperatures, a slow bath can generate considerable 
entanglement, a feature absent in the Markovian case. We mention that re- 
cently, quantum entanglement of two optical two-level systems coupled to a 
common localized environmental mode has been studied beyond the Markov 
approximation at zero temperature [9j. 
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2. Model for a single chromophore pair 

A single chromophore (index i) is modeled as a quantum two level system 
described by Pauli matrices r^ y z with energy gap Ei between ground-state 
\gi) and excited state |ej) j^]. The protein-solvent environment is formalized 
as a harmonic bath yielding the standard spin-boson Hamiltonian for each 
chromophore |5j. The Forster coupling between two chromophores is given 
by Hn = ^r( T x T x + T y T y) • Observing that the fluorescence lifetime of 
the chromophores is much larger than all other time scales, we may neglect 
the radiative decay Then, the system's total number of excitations is a 
constant of motion, and the two-chromophore system can be reduced to a 
single spin-boson model of one chromophore pair (with the Pauli matrices 

H x = ya z + ^a x + ha z X *( b l + M + > W 

K K 

where e = E\ — E 2 , and b K describe bosonic bath operators with couplings 
X K . We consider equal chromophores E\ = E% the effective basis for a 
chromophore pair is given by {| f) = \z\Q-i), I I) = | f?ie 2 ) } - 

The spectral density 1(| follows from a microscopic derivation 11]. Dif- 
ferent forms of a Debye dielectric can be assumed, but in any case, lead to an 
Ohmic form, G(ui) = 2nauje~ ul ^ ulc . The dimensionless damping constant a of 
the protein-solvent can be related to the parameters of the dielectric model 
11]. One finds for the order of magnitude of a ~ 0.01 — 0.1 0, [llj. We use an 



exponential form of the cut-off at frequency uj c . This sets the time-scale for 
the bath dynamics and is related to the reorganization energy E r ~ 2ahw c . 
If A oj c and a C 1, the bath evolves fast compared to the system and loses 
its memory quickly, rendering a Markovian approximation and the standard 
Bloch-Redfield description [12| suitable. This situation is ubiquitous in many 
physical systems, e.g., ion traps, quantum dots, and superconducting devices 



13(. It is qualitatively different for the excitons in the biomolecular environ- 
ment. Here, hw c is typically of the order of ~ 2 — 8 meV, while the Forster 
coupling constants can range from HA ~ 0.2 meV-100 meV 0,113 • Hence, 
the bath responds slower than the dynamics of the excitons evolves, and 
non-Markovian effects become dominant. Coherent oscillations in a strongly 
damped two-state system with a > 1 and A > uj c have been found using 



quantum Monte Carlo simulations [14j , Il5l | and by applying the numerical 



renormalization group 16 



3 



P(t)0 




Figure 1: Population difference P(t) for a single chromophore pair and full cross-over 
from a Markovian to a non-Markovian bath. Parameters are fceT = O.lhA, and a — 0.1. 



3. Population dynamics of a single chromophore pair 

Here, we use the quasiadiabatic path-integral (QUAPI) 0, Q, [sj to calcu- 
late the population difference P(t) = (u z ) t [lOj] in the single pair. We choose 
P(0) = 1. Figure [1] shows the results for a = 0.1 (parameters correspond to 
measured values summarized in Ref. [Ej]). P(t) decays with time in an oscil- 
latory way. The decay occurs faster for large u c while for small u c , the slow 
bath sustains more coherent oscillations. In general, for smaller uj c the spec- 
tral weight of the bath modes around the system frequency A is suppressed 
and the decohering influence is reduced, yielding prolonged coherence. In 
fact, choosing a = 0.1, u c = 0.1A (consistent with we find a coherence 
time of 1 ps which agrees well with the measured value of at least 660 fs 
[H, given the complexity of the setup. In passing, we note that we have 
compared with the Born-Markov result [lo[ for P(t) for small to c and found 
strong disagreement, as expected (not shown). 

4. Model for two coupled chromophore pairs 

Next, we address entanglement between two chromophore pairs under 
the influence of a slow bath. We consider two equal pairs described by a x /z,i, 
coupled by an interpair Forster interaction J and coupled to a harmonic 
bath. The total Hamiltonian reads 

H 2 = }] —^-Ox,i + hJ (&x,l<Tx,2 + O y ,lOy,-2) (2) 
i=l,2 
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Figure 2: Time evolution of the populations p^(t) for two coupled chromophorc pairs for 
a slow (w c = A) and a fast (uj c = 10A) bath, for k B T = Q.lhA, a = 0.1, and J = 0.1A. 



~ K K 

whose basis refers to the states { | |i) = ViQ-i), I li) = |#ie 2 ), | T2) = | e 3 ^ 4 > , | | 2 ) = 
|(7 3 e 4 )}. As before, the bath spectral density follows from a Debye dielectric 
model, again yielding the Ohmic form. The time-dependent reduced density 
matrix p(t) is computed using an adapted QUAPI scheme. Figure [2] shows 
the time-evolution of the populations P]]{t) = Pn(t) and p-\i(t) = Pii(t) 
of the four basis states for different values of u c for the initial preparation 
iVto) = (I T1I2) + I I1T2))/ a/2- After a transient oscillatory behavior, the sta- 
tionary equilibrium values are reached. The corresponding decay occurs on 
shorter times for large u c , i.e., fast baths, compared to the rather slow decay 
for small uj c . 



5. Entanglement of two chromophore pairs 

To quantify the two-pair quantum correlations, we study the entangle- 



ment along the negativity N(t) = max{0, —2( m - m (t)} [17|, Il8l ]. where Cmin(^) 
denotes the smallest eigenvalue of the partially transposed reduced density 
operator with the matrix elements Pmunu = Pnn,mv A separable state has 
iV = 0, while for a maximally entangled state, N — 1. 

Figure [3](a) shows N(t) for two values ou c = A, and u c = 50A, for the 
maximally entangled initial state | V^o) - Starting from N(0) = 1 we observe a 
decay to zero with small oscillations superimposed. For the Markovian bath 
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Figure 3: (a) Time evolution of the negativity N(t) for the cut-off lo c = A (black) and 
Lu c = 50A (red) for a — 0.01 and for J = 0.1 A (main) and J = A (inset). Moreover, 
ksT = O.lfoA. (b) Decay constant T as a function of the cut-off lo c for different values 
of a (symbols with solid lines) for fcsT = O.lfi-A, J = O.f A. Dashed line: Corresponding 
one-phonon result T = Toe^ A ^ c , where the proportionality constant Tq = 0.041 A has 
been obtained from a fit to the three data points ui c — 25 A, 30 A, 50 A. 



u c = 50A, the decay occurs faster than for the non-Markovian bath u c = A, 
indicating that entanglement survives on a longer time scale for the slow bath. 
For a larger interpair coupling J = A, the superimposed oscillations are more 
pronounced (Fig. El^a) inset) which is due to constructive interference of the 
transitions within the chromophore. For a quantitative picture, we fit an 
exponential N(t) = Nq exp(— Ft) + N\ with a constant T which contains the 
influence of the bath. Figure Mjo) shows the dependence of T on uj c . Clearly, 
T strongly decreases for small u> c , while for large uj c , the decay rate saturates 
to a constant value. The dependence of T on uo c is more pronounced for 
larger values of a. This observation indicates that entanglement could be at 
least as robust in biomolecular systems as in other macroscopic condensed- 
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Figure 4: (a) Negativity N(t) as a function of u> c for J = 0.1A (main) and J = (inset), 
for a — 0.01 and ksT = O.lftA. (b) Same as in (a), but for the strong coupling case 
a = 0.1. 



matter devices [13] which display quantum coherent behavior. We note that 
this tendency is already captured by the result of a one-phonon perturbative 
analysis, i.e., in second order perturbation theory. It would predict that for 
uo c < A no bath modes for efficient one-phonon processes are available and 
T oc G(A) oc e~ A / Wc . A fit to the data points in the region of large uj c is shown 
by the dashed line in Fig. [3](b) for a = 0.01. Clearly, when u c < A, multi- 
phonon contributions become significant. This is even more pronounced for 
larger a (not shown). 

To study the cross-over between fast and slow baths, we show N{t) for 
varying uj c in Fig. H] for the initial state \ipi) = | t i T2) • Figure H](a) shows 
the result for J = 0.1 A. The entanglement is rather quickly destroyed in the 
regime u c ^> A. On the other hand, we find a regular oscillatory decay for 
0.1 A < uj c < A. In this regime, complete entanglement disappearance and 
revivals alternate. The time scale of the entanglement oscillations is given by 
27r/J. The constructive role of a slow bath is further illustrated in the inset 
of Fig.H^a), where N(t) is shown for J = 0. In fact, in the regime uo c < A, we 
find that entanglement between the two pairs is generated by their common 
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interaction with a slow bath. Most interestingly, for u c = 0.1A, N(t) steadily 
grows even over rather long times up to tA = 500. In view of the single- 
pair results described above, this seems counterintuitive since for small u c , 
a reduced influence of the bath modes would be expected. However, in this 
regime, the bath is rather efficient in generating entanglement. This feature 
survives even for larger values of a, see Fig. E^b). The oscillatory behavior of 
the entanglement generation still occurs for J = 0.1 A, where N(t) assumes 
all values between zero and one. The bath-induced destruction happens here 
earlier due to the large a. Entanglement is also produced when J = 0, see 
inset of Fig. H^b), for 0.1A < uo c < A. Also here, N(t) can even reach the 
maximal value at intermediate times. 

The generation of entanglement can be qualitatively understood by per- 
forming a polaron-like transformation U = exp[z(cr 2j i + a z> 2)p/2] with p = 
^(b^ — b^). Setting J = and assuming that the two qubits are spatially 
close enough, the resulting Hamiltonian Hi = WH^U acquires an effective 
direct coupling i^int = ^^cfr cr z,i cr z,2 with J cS = —au c /8 ~ E r (note that 
restrictions to the applicability of U apply. Details will be published else- 
where). Then, the long- wave length bath modes are efficient in generating 
coherent coupling, and thus entanglement. Its dynamical generation occurs 
on a time scale 1/J e g, see insets of Fig. HI On the other hand, damping 
destroys coherence on a time scale T related to u c , see Fig. [3j i.e., for large 
u c , damping beats entanglement generation. 

So far, we have studied not so high temperatures, similar to the experi- 
mental conditions in Refs. [H, 0). However, in Fig. [5](a) (main) we plot N(t) 
for varying u c , for ksT = hA, for the initial state j^i)- We still find large 
entanglement oscillations at short to intermediate times, for 0.1 A < uj c < A 
despite the rather large temperature: this is an outstanding hardware feature 
that could provide a useful resource for the artificial design of controlled, 
robust, and efficient biomolecular nanostructures for quantum information 
processing 



19, 20, 21 



Furthermore, we have varied the initial preparation to the state ^2) = 
a| T1I2) + b\ I1T2) with a 2 + b 2 = 1. The inset of Fig. E](a) shows N(t) for 
varying a 2 and J = 0.1 A. \ip 2 ) is maximally entangled for a 2 = 1/2, for 
which N(t) decays monotonously with time, while away from this region the 
negativity again shows collapses and revivals. For the borders a 2 —>■ 0, 1, ^2) 
is a separable state, but entanglement is rather quickly generated with time 
before it finally dies out. Robust entanglement thus depends on the initial 
preparation and is favored by the choice of initially separable (or weakly 
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Figure 5: (a) Negativity N(t) for varying cu c ; J — 0.1A, a = O.Ol and fc^T = HA (Main). 
Inset: N(t) for different initial preparations \ip2) = a\ T1J.2) + b\ J.1T2), for J = 0.1A, a — 
0.01, fe s T = O.lftA, and w c = A. (b) N(t) for varying J; a = 0.01, fcsT = 0.1ftA,cj c = 
0.1A. J c = l/\/2 marks the border above which the initial ground state belongs to a DFS. 



entangled) states. 

Finally, we analyze the dependence on the interpair coupling J. The 
negativity N(t) is shown in Fig. [5](b) for varying J for the respective ground 
state as the initial preparation. From Eq. ([2]) it follows that a critical value 
J c = l/y/2 exists such that for J > J c , the state \ip g ) = (| T1I2) — I |i|2))/\/2 
is the two-pair groundstate, which, however, belongs to a decoherence-free 
subspace (DFS) of H2 22j- Hence, N(t) remains constantly maximal. For 
J < J c , the ground state has some weight outside of the DFS and hence 
suffers from decoherence. 

We emphasize that we have formulated our approach in quantum statis- 
tical terms. It directly involves the reduced density operator which describes 
mixed ensemble states. Thus, the entanglement dynamics reported here 
manifests itself in a statistical many-particle ensemble and, hence, would 
appropriately allow the design and implementation of robust biomolecular 
entanglement proof-of-principle experiments with current technology. 
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6. Conclusions 



In this Letter, we have shown that non-Markovian effects are vital for the 
correct description of the entanglement dynamics of biomolecular excitonic 
qubits. Recent results obtained from a reduced hierarchy equation approach 



231 ] are within the lines of our findings. While our QUAPI technique yields 



numerically exact results, the approach developed in [23J allows to include 
non-Markovian corrections only up to a certain extent. In either case, the 
results yield that the validity of the Born-Markov approximation for the 
treatment of slow protein-solvent environments in photosynthetic complexes 



24| is questionable, even in the regime of high temperature; see, e.g., our 



Fig. O 



Our results are relevant to molecular architectures 25, 261 and ultrafast 



processes |27J|. They could prove crucial in the design of artificial light har- 



vesters for robust biomolecular entanglement, with enhanced energy transfer 



rates 



281 ] for the control [21] of quantum bits. Most importantly, our pre- 
dictions could be tested with currently available experimental techniques 
a a a a. For instance, advanced ultrafast spectroscopy techniques allow 
to measure the time dependence of the elements of the two-particle reduced 
density matrix from which an entanglement measure can be derived straight- 
forwardly. Similarly, such an experiment would give direct insight into the 
interpair coupling strength which can directly be derived from the oscillation 
period of the entanglement measure. It should also be possible to prepare 
desired initial states. A further crucial topic is the behavior of this effect at 
higher temperature which can be experimentally doable. We hope that our 
results can be experimentally demonstrated in the near future. 
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